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Abstract 



> 

I/") ' We study the classical limit of quantum mechanics on graphs by intro- 

ducing a Wigner function for graphs. The classical dynamics is compared 
to the quantum dynamics obtained from the propagator. In particular we 
consider extended open graphs whose classical dynamics generate a diffu- 
sion process. The transport properties of the classical system are revealed 
in the scattering resonances and in the time evolution of the quantum 
system. 



1 Introduction 

<—| ' In this article we study quantum properties of systems which present transport 

O , behavior such as normal diffusion in the classical limit. For classical systems, 

^J ' transport phenomena has been related to dynamical quantities by the escape 

rate formalism fij|. In this formalism, the escape rate given by the leading 
Pollicott-Ruellc resonance, determines the diffusion coefficient of the system in 
the large-system limit. Since the classical time evolution is a good approxima- 
tion of the quantum evolution in the so-called semiclassical regime, one expects 
that kinetic properties such as the escape rate and the diffusion coefficient will 
emerge out of the quantum dynamics. Connections between the quantum scat- 
tering resonances and the classical diffusive behavior are also expected because, 
for open quantum systems, the quantum scattering resonances determine the 
time evolution of the wavefunction. 

The purpose of the present paper is to explore these kinetic phenomena in 
model systems known as quantum graphs. These systems have similar spectral 
statistics of energy levels as the classically chaotic systems [|[ ||. Since the 
pioneering work by Kottos and Smilansky, several studies have been devoted 
to the spectral properties of quantum graphs (j, |[ || and to their applications 
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in mesoscopic physics 0. We have studied the level spacing distribution of 
quantum graphs getting some exact results in simple cases J6|. They have also 
provided the first model with a semiclassical description of Anderson localization 
[pi . Moreover, the classical dynamics on these systems has been studied in detail 
in Ref. H where we have introduced a time-continuous classical dynamics. In 
this way, the quantum and classical dynamics on graphs - and the relationships 
between both - can be studied on the same ground as in other systems like 
billiards for instance. In this article, we go further in the connection between 
the quantum and classical dynamics by showing that the classical dynamics of 
Ref. H emerges out of the quantum dynamics introduced in Refs. [g, ||. 

With this correspondence established, and thanks to their simplicity, the 
quantum graphs turn out to be good models to study the quantum properties 
of systems which present transport properties such as diffusion in the classical 
limit. These properties have been previously studied in systems like the kicked 
rotor which has a classical dynamics (given by the standard map) presenting 
deterministic diffusion for some values of the parameters. In particular, the 
decay of the quantum staying probability, for an open version of the kicked 
rotor, has been compared to the classical decay obtained numerically from the 
simulation of trajectories of the corresponding open standard map pi]. It has 
been argued that the decay observed in the quantum staying probability is 
determined by the Pollicott-Ruelle resonances but no direct evidence has been 
reported. Here, we present results that show that this is indeed the case. The 
continuous time evolution of the wavefunction for graphs is obtained from the 
propagator which is obtained from the Fourier transform of the Green function 
on the graph. We present here a derivation that allows us to compute this 
Green function and, therefore, the propagator. On the other hand the classical 
dynamics on graphs developed in Ref. || allows us to compute the Pollicott- 
Ruelle resonances. The Pollicott-Ruelle resonances are of special importance 
because their quantum manifestation has been found in experimental measures 
of some correlation functions in microwave scattering [|lfj . 

Apart from the aforementioned time-dependent quantities we have also stud- 
ied spectral quantities like the quantum scattering resonances. The quantum 
scattering resonances have also been studied for the open kicked rotor. The 
distribution of their imaginary part has been conjectured to be related to the 
diffusion process observed for the classical system |lSq| . Here, we present nu- 
merical support for this conjecture by showing that the widths of the quantum 
resonances have the power-law distribution of Ref. |12| for some multiconnected 
diffusive graphs. 

The article is organized as follows. In Section 0, we define the quantum 
graphs and we review some of their main already known properties, namely, 
the formulation of their quantization and its exact trace formula. Section 
presents the problem of scattering on quantum graphs. In Subsection [f.2| , we 
introduce a multi-scattering expansion for the Green function on graphs. Green 
function on graphs have been considered elsewhere H, but to our knowledge 
the multi-scattering expression and the resumed closed form that we obtain are 
new results. The knowledge of the Green function allows us to obtain the prop- 



agator on graphs and therefore to have access to time-dependent phenomena 



on graphs. The propagator is introduced in Subsection 4.3. The emergence of 
the classical dynamics out of the quantum dynamics is studied in Section H. 



In Subsection 5.1, we introduce a Wigner function for graphs and we compute 
the classical limit by neglecting interference between different paths. This limit 
corresponds to the classical probability density on the graphs as we show in 



Subsection 5.2 where we also summarize the most important results about the 
classical theory of graphs |9. The quantum time evolution of open graphs is 



considered in Subsection 5.4 in which we compare the decay of the quantum 



staying probability with the classical decay of the density as obtained from the 



Pollicott-Ruellc resonances presented in Subsection 5.5. We do this comparison 
for small systems and for large systems where a diffusion process dominate the 
classical escape. In Section g, we analyze the statistical properties of the distri- 
bution of quantum scattering resonances. Here, we present a new derivation of 
the resonance density using the concept of Lagrange mean motion of an almost- 
periodic function. Examples are given and discussed in Section fij. Conclusions 
are drawn in Section H. 

2 Quantum graphs 

2.1 Definition of graphs 

Let us introduce graphs as geometrical objects where a particle moves. Graphs 
are V vertices connected by B bonds. Each bond b connects two vertices, i and 
j. We can assign an orientation to each bond and define "oriented or directed 
bonds". Here, one fixes the direction of the bond [i,j] and call b — (i,j) the 
bond oriented from i to j. The same bond but oriented from j to i is denoted 

b — (j, i). We notice that b = b. A graph with B bonds has IB directed bonds. 
The valence Vi of a vertex is the number of bonds that meet at the vertex i. 

Metric information is introduced by assigning a length If, to each bond b. In 
order to define the position of a particle on the graph, we introduce a coordinate 
Xf, on each bond b = [i,j]. We can assign either the coordinate Xuj\ or x/j t {\. 
The first one is defined such that xuj) = at i and xuj) = h at j, whereas 
x U-i) = at j and x/j,i) = h at i. Once the orientation is given, the position of 
a particle on the graph is determined by the coordinate Xf, where < Xf, < If,. 
The index b identifies the bond and the value of x\, the position on this bond. 

For some purposes, it is convenient to consider b and b as different bonds 
within the formalism. Of course, the physical quantities defined on each of 
them must satisfy some consistency relations. In particular, we should have 
that l~ b — If, and xi = If, — x&. 

We introduce here some notations that we are going to use next. For oriented 
bonds b we define the functions q(b) and p(b) which give the vertex at the origin 
and at the end of b, respectively. Thus, for the bond b — (i, j), we have q(b) = i 
and p(b) — j. These functions are well defined for graphs with multiple loops 
and bonds also. In the last case these functions take the same values for two 



or more different bonds. Note that we have the following equalities p(b) = q(b) 
and q(b) = p(b). 

3 Quantum mechanics of a particle on a graph 

On each bond b, the component ipb of the total wavefunction $ is a solution 
of the one-dimensional Schrodinger equation. This means that the dimension 
of the vector W = [ipi (#),... , ■0b(zO] T is B, but when we consider directed 
bonds as different it will be of dimension 2B, with B components containing 
redundant information [see Eq. (y) below] . We consider the time- reversible case 
(i.e., without magnetic field) 

d 2 
--T-2^b(x)=k 2 ip b {x), b=(i,j), (1) 



where k = \j2mE /% is the wavenumber and E the energy. [We use the short 
hand notation ipb(x) for ipb(xb) and is understood that x is the coordinate on 
the bond b to which the component ipb refers.] Moreover, the wavefunction must 
satisfy boundary conditions at the vertices of each bond (i and j in the previous 
equation). The solutions will have the form 

ipb(%) = ip+(b) exp(ikx) +ip-(b) exp(— ifcir) (2) 

where the boundary conditions impose restrictions on ?p+{b) and tp-(b) which 
are the amplitudes of the forward and backward moving waves on the bond b. 
If we consider oriented bonds, the wavefunction in a bond b and in the 
corresponding reverted bond b must satisfy the following consistency relation 

i/>l(x)=ip b (lb-x)- (3) 

3.1 Boundary conditions and quantization conditions 

3.1.1 Vertex matrix 

A natural boundary condition is to impose the continuity of the wavefunction 
at all the vertices together with current conservation, i.e., 

lpb(0) = <y3i 

for all the bonds b which start at the vertex i and 

4>b(h) = <Pj 
for all the bonds b which end in the vertex j. The "current conservation" reads 

^5[q{b),i}—^ b {x) =-aci(pi (4) 

LLJu x — ^U 



where J^b denotes a summation over all the IB directed bonds. The Kronecker 
delta function selects those bonds which have their origin at the vertex i and 
we impose this condition at all vertex i. The case where a, = is referred to as 
Neumann boundary condition. For nonvanishing on, Eq.(H) is the most general 
boundary condition for which the resulting Schrodinger operator is self-adjoint 
[Ol [ll|. In the case where the vertex i connects two bonds, cti plays the role of 
the intensity of a delta potential localized at the position of the vertex. In fact, 
Eq. (|j) represents the discontinuity of the derivative of the wavefunction at the 
point where the delta potential is located. Due to this analogy we call the at 
the vertex potentials. 

This boundary condition gives ||, [F|] 

J2 KVi = ° ( 5 ) 

where the sum is over all the vertices of the graph and the vertex- vertex matrix 
h is defined by 

This expression for h is valid for graphs with loops and multiple bonds connect- 
ing two vertices. Note that a loop contributes twice to each sum. In the absence 
of loops and multiple bonds it simplifies to 

f S(i,j)J2 b {6[q(b),i] cotkl b -f} , i=j 

h l3 = l i (6) 

I ~^ij sinkl b ' l r 3 

where Cij is called the connectivity matrix. Eq.(ra) has solutions only if 

deth(fc) = . 

This is the secular equation that gives the eigenenergies {fc^} of the graph. This 
equation can be solved numerically to obtain an arbitrary number of eigenen- 
ergies. Therefore, graphs are very nice systems to study the statistics of the 
spectrum fluctuations. 

3.1.2 Bond matrix 

There is however a general boundary condition, where we consider a scattering 
process at each vertex. In each vertex a (unitary) scattering matrix relates 
outgoing waves to the incoming ones. If we call a 1 the scattering matrix for the 
vertex i the condition is 

€ ut «=E^ n (i) 

b 



where the sum is over all the (non-directed) bonds that meet at i. For a ^ b, 
a l ao is the transmission amplitude for a wave that is incident at the vertex i 
from the bond b and is transmitted to the bond a. Similarly, a aa is a reflection 
amplitude. The matrix a % has dimension i/, x Vi where vi is the valence of the 
vertex i. This equation is imposed at all vertices and for all the bonds that meet 
in the vertex. 

Now we consider oriented bonds which allows us to write the previous bound- 
ary condition in the following way 

r a ut Ma)} =J2°a b a) 5{q(a),P(b)W{p(b)} (7) 

b 

where the sum is over all the 2B directed bonds and the equation is imposed 
on every directed bond a because the set of all directed bonds is equivalent to 
the set of all vertices with the (non-directed) bonds that meet at the vertex. 
As a consequence of Eq. (0) we have the relation 

i,?\p(b)]=eMikhW nt [«(&)]• (8) 

Setting the expression (||) in Eq.(0) we get 

C ut [«(«»)] = Yl T ab eMikl b W nt [«(&)] 

b 

from which follows the quantization condition 

det[l - R(fc)] = (9) 

with 

R^TD(fc) (10) 

a unitary matrix of dimension 2B where 

Dab = 8 ab e kW , with l a = l b (11) 

and 

T ab = a q ab a) S[q(a),p(b)} (12) 

Eq.(ph gives the eigenenergies {^} . Note that T ab is the transmission amplitude 
from b to a [if they are oriented such that p(b) — <?(&)]■ The reflection amplitude 
is now T aa and not T aa that vanishes due to the delta function in Eq.(12|). 

The general boundary condition (|7]) reduces to the previous one (|4) for the 
choice 



n 



S ab ) 5\p(b),l}6[q(a),l} (13) 



with @ 

uji = 2arctan ( t— ) . (14) 

\kvij 

From Eq.(pf) an exact trace formula can be obtained pi a 

d(k) = — + - V AU P cos(rfcL) (15) 

p.r 

where Ltot/i" gives the mean density of levels and the oscillating term is a 
sum over prime periodic orbits and their repetitions. A p — T a b ■ ■ ■ T za is the 
probability amplitude of the prime periodic orbit and plays the role of stability 
factor including the Maslov index. The Lyapunov coefficient per unit length of 
the orbit A p is defined by the relation \A p 



2. _ gAptp 



4 Quantum scattering on graphs 

So far we have considered bonds of finite lengths. When we attach semi-infinite 
leads to some vertices, the physical problem changes because there is now escape 
from the graph and, thus, it must be analyzed as a scattering system. The 
scattering matrix is a square matrix of a dimension that equals the number of 
open channels. For graphs these channels have a concrete meaning, they are 
the L semi-infinite leads of the graph. In this section, we shall introduce the 
scattering matrix S for graphs and we shall show that it has a multi-scattering 
expansion, closely related to the one we shall obtain for the Green function. 
Scattering on graphs was also studied by Kottos and Smilansky who showed 
that they display typical features of chaotic scattering [jl6| . 

4.1 Scattering matrix S 

On each bond and lead the Schrodinger equation allows counter propagating 
solutions. We denote by c the scattering leads and by c = {c} the set of all 
the scattering leads. Moreover, we denote by b = {b} the set of all the bond 
forming the finite part of the open graph without its leads. The Lx L scattering 
matrix S relates the incoming amplitudes to the outgoing ones as 

* out (c) = S * te (c) (16) 



We shall derive the matrix S starting from Eq.(m) which together with Eq.(|l2|) 
reads 

c ut = Ewi, n a?) 

b 

This equation is valid for every directed bond and every directed lead. We have 
dropped the explicit dependence on the vertex in equation (|17|) [see Eq.(pl)] 



because it is always understood that the incoming wave tp™ in the directed 
bond b is incident to the vertex p(b) while the outgoing wave tp% ut emanates 
from q(b). This convention is used through all the present paper. We assume 
that the scattering leads are oriented from the graph to infinity. Since the leads 
are infinite there is no scattering from the lead to the reverted lead that is 
the L x L sub-matrix Te c = 0. Neither there is transmission from a bond to 
a reverted leads nor from a lead to a bond. That is the L x 2B sub-matrix 
Tab = and the 2B x L sub-matrix Tb c = 0. Moreover, Tee = and T cc = 
due to the delta function in definition (n2|) and the selection of the orientation 
of the leads. Thus, in matrix form, Eq.(|l7|) reads 



*out(c) " 












" *in(6) 


*out(c) 


= 


1 cc 


T cb 




*in(c) 


*out(b) 




Tbe 


T bb 




*in(b) 



where T c e is a L x L matrix whose elements represent the direct lead-to-lead 
transmission or reflection amplitudes. The matrix T cb is a L x 2B matrix whose 
elements represent the bond-to-lead transmission and, similarly, the matrix T b e 
is a 2B x L matrix whose elements represent the lead-to-bond transmission. 
Finally, the matrix T bb is the 2B x 2B matrix that we have simply called T 
for the closed graphs and represents bond-to-bond transmission or reflection. 
Note however that, for a vertex with an attached scattering lead, the bond- 
to-bond probability amplitude is different from the one for the graph without 
the scattering lead because the valence of the vertex has changed from Vi to 
V{ + ^{lcads attached to i}. The previous matrix equation can be rewritten as 



*out(c) = T c e* in (c)+T cb * in (b) 
tfout(b) = T b e* in (c)+T bb * in (b) 



(18) 
(19) 



and now we use the relation given in Eq.(0). With the convention introduced 
after Eq. ( J17| ) , we can write Eq. (@) in a matrix form 

% n (b) = D(fc)* out (&) (20) 

with the diagonal (2B x IB) matrix D defined in Eq.©. Thus from Eq.© 
and Eq. (po|) we get 

* in (b) = (D- 1 - T bb )- 1 T b e* in (c) . 

Replacing this last equation in Eq.dlS) we obtain 

tfout(c) = {T c e +T cb [l - R(fc)]- 1 DT b e}* i „(c) 

where we wrote R(/c) = D(fc)T bb . That is the outgoing waves on the leads are 
determined by the incoming waves on the leads. This gives the desired scattering 
matrix 



S = T c6 + T cb [l-R(A : )]- 1 DT b e 



(21) 



which appears in ( pf ) once we identify each lead and its reverse as the same 
physical lead. The multiple scattering expansion is obtained from 

oo 

n-R(k)]- i = Y,m n - 

ra=0 

In a similar way as was done for the trace formula we can get that 
S CC ' = [T C c] cc> + ^2 A p exp[ikl p (c',c)} 

with V c '^c the set of trajectories which goes from c' to c. As usual A p is the 
amplitude of the path and l p (d , c) the length of the path given by the sum of 
the traversed bond lengths. 

4.2 Green function for graphs 

Following Balian and Bloch [17) we seek for a multi-scattering expansion of the 
Green function. The Green function represents the wavefunction in the presence 
of a point source. We shall identify the bond where the point source is with the 
bond 1. In this bond, the point source is located at X\ — x 1 . Since we shall 
work with directed bonds the source also appears at the point Xf = h — %' on 
the bond 1. The Green function satisfies the following equation 

d 2 2m 

— xGi(x,x') + k Gx(x,x') = —rS(x - x') , < x,x' < h 
ax 2 - h 



d 2 2m 

—^Gi(x,x') + k 2 Gi(x,x') = —5(x- h + x') , < x < Zj = h 
ax h. 



d 2 
-^G b (x,x') + k 2 G b (x,x / )=0, V6^{1,1} 

Note that x' is fixed to belong to the bond 1. The Green function satisfies the 
same boundary condition as the wavefunction. 

In the spirit of a multi-scattering solution we assume that 

Gb(x, x 1 ) = G (x, x')5 b i + G (x, h - x')S- bl + 

Go(x, 0)n q (b)(x') + G (x, k)Hp(b)(x') (22) 

where Gq(x,x') — ^p£ e 2ik — is the free Green function and represents an out- 
going wave from x' . The justification of this ansatz is the following. In the 
bond 1, the wavefunction consists in the superposition of the wave emanating 
from the source at x' plus waves that arrive from the borders of the bond. On 
the other bonds (i.e. not 1 or 1) only these waves are present. Since we are 



dealing with directed bonds we call the vertices at the border of b, q(b) and p(b) 
as usual. 

We have to impose to the Green function the consistency condition 

G b (x,x') = G b (l b -x,x r ). (23) 

The observation that 

G (x,h - x) = Go(h — x,x') 
G o (h-x,0) = G (x,l b ) 
Go(lb — x,lb) = G (x,0) 



leads us to conclude from Eq. (|23|) that 

^q(b)( x ') = Vp{b){x') 

Mp(6)( X ') = M?(6)(*0 

Now, we impose the boundary condition which is given by Eq.dlj). With this 
aim we have to identify incoming and outgoing components. The amplitude of 
the incoming wavefunction from the bond b is 

C = G (h, x')6 bl + G (h, h - x')S bi + G (l b , 0)n qib) (x') 

where the first two terms represent the incoming amplitudes of a wave emanating 
from x' if we are on the bonds 1 or 1. The third term is the incoming amplitude 
of the wave transmitted to the origin of b. For the outgoing amplitude on the 
bond b we have 

r b ut = G (0,0) N{b) (x') 



Thus the boundary condition Eq.(17) gives 
2M?(o)(*') =^2T ab [G (lb,x')8 bl +G Q (lb,h ~x')8 b i + G Q (l b ,0)n q ( b )(x')} 



ikh 

b 



where we used that G (0,0) = -^. We define g{x,x') = l -^-G (x,x'). If we 
agree in denoting by T g (l,0) the matrix whose ab elements are T ab g(lb 7 0) and 
by T g (l,x') the matrix whose ab elements are T ab g{l b) x') we can rewrite the 
previous system of equations in the more convenient vectorial form 

H q (x') = T g (l, x') ■ ei + T g (l, h - x') ■ ej + T g (l, 0) • fi q (x') 

where ei (ej) is the 2_B-dimensional vector of components [ei]j = S b i ([ej]b = 
5 b i). The solution to these equations can be written as 

fi q (x') = [I - J g (l, 0)}- 1 [T g (l, x') ■ ei + T g (l, h - x') ■ e t ] . (24) 



10 



Replacing in ( p2| ) with fj, p ( b )(x') — M q (h)( x ')> we nave obtained the Green func- 
tion for graphs. We note that T g (l, 0) = R(fc) and thus the Green function 
has poles at the resonances. The multiple scattering form follows from the 
well-known expansion (I — R) _1 = I + R + R 2 + R 3 + • ■ • . Remembering that 
H p t a )( x ') = ^q(a)( x ')i we have that the a component of the Green function is 

G a {x, x') = G {x, x')S al + G (x, 0) \T al g{h,x') + T ai g(0, x')} (25) 

+ G (x, l a ) [T &1 g(h,x') + T ai g(0, x')] 

+ G (x, 0) Y T ab g(l b , 0) [T bl g(h, x') + T bi g(0, x')] 

b 

+ G {x, l a ) Y Tabg(k, 0) [T bl g(h , x') + T bi g(0, x')] 

b 

+ G {x, 0) Y T ab9(h, 0)T w g(l b ,,0) [T vi g{h,x') + T vl g{Q, x')} 

b,b' 

+ G {x, l a ) Y T &bg{k, 0)T bb ,g{l b ,,0) [T bll g{h,x') + T b , i9 (0, x')} + ■■■ 



b.b' 



Noticing that g(x,x') — e lk \ x x ' we can write the previous result in the more 
handy notation 

G a (x, x') = ^±Y A P e ikl »^ (26) 

M 

where A p is the probability amplitude of the path p which connects the initial 
point x' on the bond 1 to the final point x on the bond a. If the path p is 
composed by the n bonds 1& 2 • • • b n -ia then 

A p = T a b rl _ 1 T bri _ 1 i> n _ 2 ■ ■ -T b2 i 

The fact that Gq(x,x') and g(x,x') depend on the modulus of the differences 
\x — x'\ implies that, in (p5[), we are always adding lengths. Thus l p (x, x') is the 
total length of the path s which connects x 1 to x. 

The expression (E6|) is like a path-integral representation of the Green func- 
tion: We add the probability amplitudes of all the paths connecting x' to x in 
order to get the Green function. 

4.3 Propagator for graphs 

Green functions and propagator are related by Fourier or Laplace transforms: 

G(+\x,x';E) = lim — / dtexp(-et)exp ( ~Et ]K a (x,x';t) 
e ^o+ in J \h J 

Ka(x,x , ;t) = ^- f dEexp(-~Et\G a (x,x';E) (27) 



11 



where the contours C+ goes from KeE = +00 to RcE = —00 with a positive 
imaginary part, while C_ goes from Re E = —00 to Re E = +00 with a negative 
imaginary part. 

Using expression (p6[) for the Green function, we get from Eq.(p7|) that 



^*<<) = V^£^"^ (28) 

v M 

This expression shows that the propagator is the sum over all the paths {p} 
that join x' to x in a fixed time t. Each term is composed of a free propagator 
weighted by the probability amplitude of the given path. This result could have 
been guessed from the general principles of quantum mechanics, i.e., if there are 
many ways to obtain a given result then the probability amplitude is the sum 
of the probability amplitudes of the different ways of obtaining the result. 

The closed form of the Green function [given by Eqs.(|22j) & (p4[)1 and the fast 
Fourier transform allow us to obtain numerically the propagator as a function 
of the time t and, therefore, the time evolution of a wave packet. We shall 
develop this possibility after analyzing the classical limit of quantum mechanics 
on graphs in the following section. 

5 Emerging classical dynamics on quantum graphs 

The emergence of the classical dynamics out of the quantum dynamics can be 
studied by introducing the concept of Wigner function. Such a function should 
tend to the classical probability density in the classical limit. 

5.1 Wigner functions on graphs 

The so-called Wigner function was introduced by Wigner in order to study 
systems with a potential extending over an infinite physical space. For graphs, 
we cannot use the same definition since each bond is either finite or semi-infinite. 
We define a Wigner function for a graph in the following way. On each (non- 
directed) bond a the Wigner function is given by 

fa{x,p) = — dye^Mx-y/2Wa(x + y/2) for 0<x<^ (29) 

2Trh J_ 2x 2 



and 



fa(x, P ) = -^- / X dye** M* ~ v/2)<(s + y/2) for k < x < \ a 

(30) 



In this way, the argument of the wavefunctions always remains in the interval 
(0, Ig) corresponding to the bond a. We notice that f a (x = 0,p) = f a (x — 
la,p) =0 with of the definitions (p9|) and (|30|). 



12 



The Wigner function is a "representation" of the wavefunction in phase space 
and it is essential to have a unique correspondence between the Wigner function 
and the wavefunction. In order to show that this is the case with our definition, 
we multiply Eq.(E3) by e l ~s~ and we integrate with respect to p to get 

+00 

dpe^f a (x,p) = M* + 2//2)€(z - y/2) (31) 

-00 

If we set y = we obtain the probability density on the bond a 

\M*)\ 2 = / dpf a (x,p) (32) 



On the other hand, if we set x = y/2 in Eq.(31) we get 



dpe t -f a (y/2,p) = MyWa(0) 

which is valid for < x — y/2 < y, i.e. < y < l a , and we recovered the 
wavefunction (or its conjugate) on all the bonds except for a constant factor 
that is fixed by the boundary condition and normalization. We could also have 
proceeded in a similar way with Eq.(Bfl). We notice that this result shows 
that there is redundant information in the Wigner function. The conclusion is 
that the wavefunction on all the bonds is encoded in all the Wigner functions 
{/„(z,p)}f =1 . 

Since a wavefunction can always be written in terms of the propagator (which 
is also solution of the Schrodinger equation and represents the evolution from 
a delta localized initial state) we will compute the Wigner function for the 
propagator. Other cases are obtained by convenient averages over the initial 
conditions. Thus we need to compute K(x — y/2, x'\ t)K*(x + y/2, x'; t). With 
this purpose, we use the expression ( |28|) which expresses the propagator as a 
sum over paths (in this section we use the letter s to refer to paths in order to 
avoid confusion with the momentum p) : 

jy, ^^ m [i s (x- H /2,x') 2 -! s ,(x + a/2,x') 2 l 

K a (x - y/2, x'; t)K* a (x + y/2, x'; t) = — £ A s A s ,e l ^ 



With l s (x, x') the length of the trajectory s that joins x' to x. 

We note that for a path that starts at x' on a bond 6 and ends at x on the 
bond b the lengths can be 

• l s (x, x') — (Jb — x') + l s + x, if the path goes from x' to the end of bo then 
eventually traverses other bonds adding a distance l s and arrives at the 
position x of the bond b via its origin. 

• l s {x, x') = x' + l s + x, if the path goes from x' to the origin of bo then 
eventually traverses other bonds adding a distance l s and arrives at the 
position x of the bond b via its origin. 



13 



• l s (x, x') = (lt, — x') + l s + (lb — x), if the path goes from x' to the end of 
60 then eventually traverses other bonds adding a distance l s and arrives 
at the position x of the bond b via the end of b. 

• l s (x, x') = x' + l s + (lb — x), if the path goes from x' to the origin of 60 
then eventually traverses other bonds adding a distance l s and arrives at 
the position x of the bond b via its end. 

We evaluate now the difference l s (x — y/2, x') — l s (x+y/2, x') for equal paths. 
We obtain 

l s (x - y/2, x') - l s (x + y/2, x') = -y 

l s (x - y/2, x') - l s (x + y/2, x') = +y 

The first result holds for trajectories that arrives at the final bond via its origin 
and the second result for trajectories that arrives via the end. Now we compute 
the sum which is in both cases 

l s (x - y/2, x 1 ) + l s (x + y/2, x') = 2 l s (x, x') 

Using the identity a 2 — b 2 = (a + b) (a — b) we have the results 

l 2 s (x - y/2, x 1 ) - l 2 (x + y/2, x') = -2 y l s (x, x 1 ) 

l 2 s (x - y/2, x 1 ) - l 2 s (x + y/2, x') - +2 y l s (x, x 1 ) 

for paths that arrives through the origin or the end of the bond b respectively. 
We took care of equal paths because we want to separate their contribution to 
the Wigner function which we call the diagonal term 



[K a {x - y/2, x'; t)K* a (x + y/2, x'; t)] di 



i;iK 



£ <e --^^ +£<e+-^ 



2ttM 

K si s 2 ) 

where S\ is the set of paths that arrive at the bond a through the vertex at the 
origin of the bond a and S2 are those paths that arrive through the vertex at 
the end. 

For the non-diagonal term the differences in the exponent are always of the 
form 

l 2 s (x - y/2, x 1 ) - I 2 , (x + y/2, x') = [l s (x, x') - l s , (x, x')][l s (x, x') + l s > (x, x') ± y] 



or 

l 2 (x - y/2,x') - l 2 ,(x + y/2,x') ^ [l s (x,x') - l s ,(x,x') ±y}[l s (x,x') + l s ,(x,x')} 
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hence we get non-diagonal terms of the form 



K -a\ x + y/ Z > x ; tjjnon-diag = 

(v-^i in \l = ( x , x ) —li(x,x~) m 


psCw.wOi^/ (»»»')] 











The calculation of the Wigner function (P9|)-(p0|), requires the evaluation of 
integrals of the form 

f +x ° ■, + ™« sinfi(p±m , s 

/ e l{ - p±a) %dy = 2Ti r ^f- — '- 33 

J-xo P±^ 

where f2 = m 3 ^ ,:c ' for the diagonal term and O = — — — ; — ^-^- — '■ — - for the 
non-diagonal terms. In the classical limit we have that 

sin^(p±0) , 

lim &¥— — - = irS(p ± n) 34 

from which we obtain the Wigner function in the limit % — ► 0. In the classical 
limit, the phase variations of the non-diagonal terms are so wild that the total 
sum is zero due to destructive interferences. We have thus the result that, in 
the classical limit % — » 0, the Wigner function defined for graphs becomes 



t l 4-\ m 



27rfii 



Y J Al 1 5\p + ml Sx {x)/t]+Y J A 2 82 6^-ml S2 {x)/t}\ (35) 

.81 S 2 ) 

This limit corresponds to the motion of the classical density in the phase space 
of the corresponding classical system. In the next section, we shall establish 
that, indeed, the classical dynamics evolves the probability density in phase 
space according to Eq.(|3F 



5.2 Classical dynamics on graphs 

We have computed the classical limit for the Wigner function on graphs [see 
Eq.(p5[)1, which should be solution of the classical "Liouville" equation. In this 
section, we shall summarize the main result obtained in Ref. H, where we 
studied in detail the classical dynamics on graphs, and we shall show that the 
density ( p5| ) is the solution of the classical equation. Therefore, the classical 
dynamics which we discuss here is the classical limit of the quantum dynamics 
on graphs as obtained from the classical limit of the Wigner functions. 

On a graph, a particle moves freely as long as it stays on a bond. At 
the vertices, we have to introduce transition probabilities Pw = \Tbb'\ ■ This 
choice is dictated by the quantum-classical correspondence as we shall see in 
this section. The dynamics is expressed by the following master equation (we 
consider the notation Xb = [b, x] ) : 

p([b,x],t)=Y, p »»f , (i b '> x '^ t - x+lb '~ x ' ) ( 36 ) 
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The time delay correspond to the time that take to arrive from x' in the bond 
b' to x in the bond b. 

The density p ([b, x] , t) defined on each directed bond is a density defined in 
a constant energy surface of phase space. In fact, the conservation of energy 
fixes the modulus of the momentum and, therefore, the points of the constant 
energy surface are given by the position on the bond and the direction on the 
bond, that is by position on directed bonds. 

The properties of these classical dynamics for open and closed systems are 
described in Ref. M where we also show that it can be understood as a random 
suspended flow. 

To establish the connection with the classical limit of the Wigner function 
we iterate the master equation (36). In Ref. g we have shown that iterating 
the master equation allows us to obtain the density at the current time t in 
terms of the density at the initial to — 0, which gives an explicit form for the 
Frobenius-Perron operator 



p([b,x],t) =^J ^2 p bb'Pb'b" ■■■P b (n-i )b( 

n (,'6"...f,(n) 



.,,( 



,(«) (n) 
, X 



•0 



(37) 



with 



x {n) = x 



vt + ^2 lb(i} 



s=l 



Accordingly, p ([&, x] , t) is given by a sum over the initial conditions 



o , x 



and over all the paths that connect 



, X 



to [6, x] in a time t. Each given 



path contributes to this sum by its probability multiplied by the probability 



o , x 



.0 



density p ( 

If the initial distribution is concentrated on a point, i.e., 

p( b { "\x n) ,0J =S bMb J(x* -x (n) ) 

Eq.(p7|) can be expressed as [we use the property S(ax) — 5(x)/a 

p([b,x],t) = -^2 ^2 Pbb'Pb'b" •••■PbC«-i)b(») 



n 6'fe". ..&(«) 



(38) 



JbM b 



(39) 



5.3 Connection with the classical limit of Wigner func- 
tions 

The probability density at a given time t is provided by the classical Frobenius- 
Perron operator (|37J). For the particular initial condition (BS) this leads to 



1G 



Eq.(B9h. Remembering that the probability of a path s was written as A 2 S (we 
use again the letter s to denote a path since p is used for the momentum) and 
noticing that the sum in Eq.(39) is a sum over all the paths connecting &* to b 
we can rewrite Eq.(B9) as 



p([b,x] ,t) 



1 



E 



Aid 



s(6»->6) 

and by the property of the delta function 



l s (x,X») 



(40) 



Is 

v 



t 



L s \X 7 X* J 



t 



i s \X^ X* ) 



t 



~ T 



1 1 lib g\tij . ib ]|t 



so that Eq.(|4(j) is equivalent to 



p([b,x],t) = - J2 A " s 

s(b»-*b) 



I I Lb c \ <b ■ *b ^ 



(41) 



where p — mv. Eq.(ETl) is (up to the normalization factor — ^t) the classical 
limit of the Wigner function as obtained in Eq. ( p5| ) . The Wigner function ( p5| ) 
also contains a term with — p because in Subsection 5.1 we defined the Wigner 



function for non-directed bonds although we deal with directed bonds in the 
present section. Eq. ( [4l|) is the probability density of being in the oriented bond 
b with momentum p. In the reverted bond b we have also a probability density 
which is obtained from the density in b by revering the sign of p. Therefore, the 
probability density of being in the non-directed bond b is 



m 
T 



E 



K* 



7Tilg2 \X j X^ J 

t 



in 

T 



E 

«i (6,— b) 



Aij 



ml Sx (x,x*) 



t 



with s\ and S2 the set of paths defined for the Wigner function in Subsection 
5.1 . The comparison with Eq.(|35|) shows that the quantum time evolution of the 
Wigner function corresponds to the classical time evolution of the probability 
density given by the classical Frobenius-Perron operator (|37]) in the classical 
limit: 



f b (x,p;t) ~ 2-*p([6,aj],t) 



for 



Accordingly, the classical dynamics introduced in Ref. 



(42) 



and summarized 



in Subsection 5.2 is the classical limit of the quantum dynamics on graphs of 
Section |[ 



5.4 Quantum time evolution of staying probabilities 

The preceding results shows that the classical dynamics emerges out of the 
quantum dynamics of quantities such as averages or staying probabilities which 
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can be denned in terms of Wigner functions. Since the Wigner functions have 
a Liouvillian time evolution according to the Frobenius-Perron operator in the 
classical limit Ti — > we should expect an early decay given in terms of the 
Pollicott-Rucllc resonances which are the generalized eigenvalues of the Liouvil- 
lian operator. For graphs, the Pollicott-Ruelle resonances have been described 
in Ref. ||. Our purpose is here to show that, indeed, the Pollicott-Ruelle res- 
onances control the early quantum decay of the staying probability in a finite 
part of an open graph. 

The quantum time evolution of the wave function is obtained from the prop- 
agator as 

tfr(t) = K(i) V(0) 



where ip(0) represent the initial wave packet. As we said in Subsection |L3 
the propagator can be obtained from the Green function for graphs by Fourier 
transformation. Therefore we can compute ^ b (x,t). 

The quantum staying probability (or survival probability) P(t) is defined as 
the probability of remaining in the bounded part of an open graph at time t. 
Since the particles that escape cannot return to the graph, there is no recurrence 
and the staying probability is equal to the probability of have been in the graph 
until the time t: 

p(t) = V [ b dx\Mx,t)\ 2 
b J ° 

To compute P(t) we can proceed as follows. First, we consider some initial wave 
packet ipi(y,0) on the bond 1 and with mean value E. The Green function 
G b (x,y, E) (where the index b refer to the coordinate x on the bond b and y 
is a coordinate on the bond 1 ) is computed from Eqs.(E2n and (E4). Then the 
propagator is obtained from the Fourier transform 

K b (x,y,t) = J G b (x,y,E)exp (~^ Et ) dE 

and the wavefunction is given by 



i>b(x,t)= dyK b (x,y,t)ipi(y,0) 

Jo 

P(t) is the quantum analog of the density distribution integrated over the 
graph, that is the classical staying probability. For a classically chaotic graph we 
know that the classical staying probability decay exponentially with a decay rate 
given by the leading Pollicott-Ruelle resonance (see Ref. 0). From semiclassical 
arguments, the quantum staying probability P(t) should follow the classical 
decay for short times. Indeed, using Eq.(p2), the quantum staying probability 
can be expressed in terms of the Wigner function (J2J3) according to 

P W=E/ dx / dpf b (x,p;t) (43) 

i JO J-oo 
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In the classical limit K — ► 0, Eq. (M2) implies that the quantum staying proba- 
bility evolves as 

rh r + co i 

P(t)^J2 dx dp—p([b,xlt) (44) 

For an energy distribution well localized around the mean energy E of the initial 
wave packet, we can suppose that the classical evolution takes place essentially 
on the energy shell of energy E. 

If we denote by P* the classical Frobenius-Perron operator and by po the 
initial probability density corresponding to the initial wavepacket, the quantum 
statying probability can be written as 

P(t) ~ (A|P'po> , for h -» (45) 

where we have introduced the observable A defined by 

A[b,x] = -^- for fee b 

znn 

A[c,x] = for c e c 

which is the indicator function of the bounded part of the open graph. Using 
the spectral decomposition of the Frobenius-Perron operator described in Ref. 
[pl , the quantum staying probability has thus the following early decay 

P(t)~£(^)^'(V^ ), for h^O (46) 

j 

in terms of the left- and right eigenvectors of the Frobenius-Perron operator: 

pty. = e s 3 t Vj and ptiy. = e a*ty. (47) 

(see Ref. (9J). Since the leading Pollicott-Ruelle resonance is the classical escape 
rate for an open graph, we can conclude that the quantum staying probability 
will have an exponential early decay according to 

P{t) ~ exp [-7d(v)t] (48) 



in terms of the classical escape rate 7ci(u) = ^7ci(^ = 1) where v — y / 2E/m is 
the velocity of the classical particle at the mean energy E — mv 2 /2. 

5.5 The classical zeta function and the Pollicott-Ruelle 
resonances 

We have shown in Ref. M that he Pollicott-Ruelle resonances of a classical 
particle moving with velocity v on graph can be computed as the complex zeros 
{sj} of the classical Selberg-Smale zeta function of the graph given by 

Z cl (s) = det [I - Q(s)] = . 
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in terms of the matrix 

i b i 

Qbb'(s) = Pwe~~ s ~ 

where Pw = |Tbb'| 2 are the transition probabilities. This classical zeta function 
can be rewritten as a product over all the prime periodic orbits on the graph as 



Z c i(s) = n p [l - e-^+^A (49) 

The zeros of the classical zeta function for a scattering system are located 
in the half-plane Re Sj < and there is a gap empty of resonances below the 
axis Re Sj = 0. This gap is determined by the classical escape rate which is the 
leading Pollicott-Ruelle resonance: so = — 7cl- 

5.6 Emerging diffusion in spatially extended graphs 

On spatially extended graphs, the classical motion becomes diffusive. We showed 
in Ref. [pi that the diffusion coefficient for a periodic graph can be com- 
puted from the Pollicott-Ruelle resonances of the extended system by introduc- 
ing a classical wavenumber q associated with the classical probability density. 
Therefore, the leading Pollicott-Ruelle resonance acquires a dependence on this 
wavenumber according to 

So (q) = -Dq 2 + O(q 4 ) (50) 

where D is the diffusion coefficient. 

If the spatially extended periodic chain is truncated to keep only N unit cells 
and semi-infinite leads are attached to the ends, we have furthermore shown 
in Ref. || that the classical escape rate depends on the diffusion coefficient 
according to 

7d(A0 « D ^ (51) 

in the limit N — ► oo. 

The previous results on the classical-quantum correspondence show that 
this diffusive behavior is expected in the early time evolution of the quantum 
staying probability for such spatially extended open graphs. This result will be 
illustrated in Section [?]. 

6 The quantum scattering resonances 

6.1 Scattering resonances 

The scattering resonances are given by the poles of the scattering matrix S in 
the complex plane of the quantum wavenumber k. These poles are the complex 
zeros of 

Z{k) = det[l - R(fc)] 
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This function can be expressed as a product over periodic orbits, using the 
identity lndet(l - R) = trln(l - R) and the series ln(l - R) = - J2 n >i R "/ n - 0ne 
gets 



Z ( fc ) = IT f 1 - e-^^ /2 e l(M " +7r ^ /2) j (52) 



where l p is the length of the prime periodic orbit p, A its Lyapunov exponent 
per unit length, and \i v its Maslov index. This formal expression is equally valid 
for open and closed graphs and thus their zeros give the eigenenergies in the 
first case (zeros in the real k axis) and the quantum scattering resonances in the 
second case, but, as for the trace formula, the product over primitive periodic 
orbits does not converge in the real axis. 

A remark is here in order about the difference between the quantum scat- 
tering resonances and the classical Pollicott-Ruelle resonances. The quantum 
scattering resonances control the decay of the quantum wavefunction and are 
defined either at complex energies E n or at complex wavenumbers or momenta 
k n . In contrast, the Pollicott-Ruelle resonances control the decay of the classi- 
cal probability density which is as the square of the modulus of the quantum 
wavefunction. Accordingly, the Pollicott-Ruelle resonances are related to the 
complex Bohr frequencies tv mn = (E m — E n )/h and have the unit of the inverse 
of a time. 

In this and the following section, we shall consider units where h = 1 
and 2m — 1, so that the quantum wavenumber k is related to the energy 
by E = k 2 . In these units, the width of a quantum scattering resonance is 
T„ = — 4Refc„Imfc„ and the velocity of a resonance is v n — 2Refc n , so that 
T„ = -2v„Imfc n . 

6.2 Topological pressure and the gap for quantum scat- 
tering resonances 

The chaotic properties of the classical dynamics can be characterized by quan- 
tities such as the topological entropy, the Kolmogorov-Sinai entropy, the mean 
Lyapunov exponent, or the partial Hausdorff dimension djj in the case of open 
systems. All these quantities can be derived from the so-called "topological 
pressure" per unit time P(/3; v) which we analyzed in Ref. || for graphs. 

Beside these important properties, the topological pressure also provides 
information for the quantum scattering problem. In fact, the quantum zeta 
function has a structure very similar to a Ruelle zeta function with some expo- 
nent (3 = 1/2 (see Refs. || [jl|). As a consequence, the quantum zeta function 
is known to be holomorphic for Im k > P(j3) where P{(3) is the so-called topo- 
logical pressure per unit length. Therefore the poles of the zeta function are 
located in the half plane 
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The following result can thus be deduced M, |T|, ||, |lj || : 

If P (h) < or equivalently if < rfn < 17 ^ e lifetimes {r„} are smaller 
than a maximum quantum lifetime r q and i/iere is a gap in the resonance spec- 
trum. 

If P (i) > or equivalently if \ < d^ < 1, i/ie lifetimes may be arbitrarily 
long. 

In the first case, the partial Hausdorff dimension is small (0 < g?h < ^) and 
we can talk about a filamentary set of trapped trajectories. In the second case 
with i < dpi < 1, the set of trapped trajectories is bulky. Hence, the result 
shows that a gap appears in the distribution of quantum scattering resonances 
in the case of a filamentary set of trapped trajectories. The gap is determined 
by the topological pressure at (3 — 1/2 which is the exponent corresponding to 
quantum mechanics, as opposed to the exponent (3 — 1 which corresponds to 
classical mechanics |lj| |2^] . 

The properties of the topological pressure yield the following important in- 
equality between the quantum lifetime r q and the classical lifetime t c \ = l/7 c i: 

- = -2 P (lv) < -P(l;v) = - 

where the equality stands only for a set of trapped trajectories which reduces 
to a single periodic orbit J18, nil . Accordingly, the quantum lifetime equals the 
classical lifetime only for a periodic set of trapped trajectories. On the other 
hand, the quantum lifetime is longer than the classical lifetime for a chaotic set of 
trapped trajectories. This result - which has previously been proved for billiards 
and general Hamiltonian systems in the semiclassical limit |l|, [19, EG, EU E2[ - 
thus also extends to open graphs. 

We note that it is thanks to the time continuous classical dynamics that we 
can compute the above estimate of the gap of resonances. Examples of this will 
be considered in Section 0. 

6.3 Mean motion and the density of resonances 

Beside the possibility of a gap in the distribution of the scattering resonances, we 
also want to obtain the distribution of the imaginary parts which give the widths 
of the scattering resonances. For this purpose, we first determine the mean 
density of resonances. This can be obtained analytically for a general graph 
and was done by Kottos and Smilansky who obtained a trace formula for the 
resonance density. Here, we proceed in a different way. The zeta function for a A;- 
independent matrix T (e.g. with Neumann boundary conditions), is an almost- 
periodic function of k and several results are known about their properties, in 
particular: 

The mean density of resonances H(yi,y2) (or zeros of the zeta function) in 
a strip yi < Im k < j/2 of the complex plane k = x — iy is determined from 
the number of resonances N(xi, X2, j/i, 2/2) in the rectangle (x±, X2, j/i, J/2) by the 
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following relation fl27| 

H(yi,yi) = hm = — [M(y 2 ) ~ M( yi )} 

|x 2 — a:i|— >oo X2 — X\ Z7T 

where M(y) is the mean motion of the function of x 

f y (x) = Z(k = x-iy). (53) 

The function 

h(y) = H(Q, y) = i- [M(y) - M (0+)] (54) 

gives the density of resonances with Im k < y. The total density of resonances 
is therefore given by h(pd). We shall compute this number using the general 
properties of the function Z(k). 

6.3.1 Mean motion 

An almost-periodic function / : t E R — » C 

f{t) = r(t)e**W 
with real r(t) and <p(t), has a mean motion M 

M = hm M 

if the limit exists. 

The problem of computing the mean motion has a long history and is a 
difficult problem. It was posed by Lagrange in 1781 and there are still only a 
few general results. For an almost-periodic function formed by three frequencies 
there is a explicit formula given by Bohl. Weyl proved the existence of mean 
motion for functions with a finite number of incommensurate frequencies and 
also gave a formula to compute it. 

The simplest result holds for the so-called Lagrangian case considered by 
Lagrange in his original work. We quote the result because it is important for 
what follows. 

Consider a function 

/(*) = a e Mo * + a ie lbJlt + ■■■ + a„e^"'. (55) 

If 

|oo| > Kl + M H h |o n | 

then / has a mean motion M which is M = ujo. We shall show that the density 
of resonances is determined by the mean motion of the function in Eq.(B3h in 
the Lagrangian case. 
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6.3.2 The density of resonances h(oo) 

Let us consider the expansion of the determinant involved in the zeta function 
[see Eq.(g)]. If 2B is the dimension of the square matrix R, then 

2B 

det(AI-R) = ^m 2B _,A' 
1=0 



where mo = l,mi = tr(— R) ... and m 2 B = det(— R). The secular equation is 



Sz=o rUl ~ 0- ^he general term 



l<i 1 <i 2 <---<i P <2B v 

where (— R) . ? I is the principal minor of order p obtained by elim- 

y 1 1 • • • lp J 

mating the n — p rows and columns of (— R) different from i\ ■ ■ -i p . This coef- 
ficient m p is the homogeneous symmetric polynomial of degree p which can be 
constructed from the 2B eigenvalues of R, therefore m p is of the formP] 

aj p e ^ jeJ " , 
{.'p\ 

where J p stands for a set of p different integers in the interval [1, 2B) and { J p } 
the set of elements J p . The important point to notice is that in m p there are p 
lengths. 

It is clear that the only term which involves all the lengths of the graph in 
the expansion 

f v (x) = det [I - R(x - iy)] = 1 + mi + \-m p -\ h m 2 B-i + m2B (56) 



m 2B = det R = e 2iLt ° tk det T = e 2iL ^x e 2L MV det T ^ 



IS 



For y — > oo, rri2B is exponentially larger than the remaining terms in Eq.(p6J) 
and the next leading term 77125-1, which can be estimated as 

m 2B -i ~ e* 2 ^-^^ (58) 

with A m i„ the minimum eigenvalue of T, is exponentially larger than m p with 
p = 1, ..., 2B — 2. Therefore we have 

\m 2B \ > \m 2 B-i\ ~\ h |1| for y — > 00 



1 Remember that R(fc) = TD(fc) according to Eq.(rLCJ). 
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Thus, the function f y (x) is in the Lagrangian case. The term rri2B corre- 
sponds to the term aoe luJot (here t is replaced by x) of the expansion Eq.(p5|) 
and the mean motion is given by the frequency of rri2B given by Eq.rt57f). Ac- 
cordingly, we find that AJ(oo) = 2L t ot- 

In order to compute the density of resonances h(oo) [see Eq.(|54|)] we have 
to evaluate M(0 + ). Since the upper half of the complex plane is empty of 
resonances we have that M(0 + ) = and, thus, from Eq.(p3) we get the density 
of resonances 

fc(oo) = — (59) 

This result prevails as long as detT ^ 0. In the opposite case the density is 
given by 

■ i \ *i'tot — 'min iai\\ 

Hoo) = — (60) 

if the corresponding constant factor does not vanish. 

From this argument, it is clear that, once the function f y belongs to the 
Lagrangian case (that is, when y = y ma x), no further resonance appears below 
2/max- This allows us to estimate how deep in the complex plane lies the shortest 
living resonances. With this aim, we consider the largest terms of the expansion 
"56b, i.e., 



f y (x) = m 2 B + m 2 B-i H 

with m 2B = e 2lLt ° tfc detT and m 2B -i ~ e i ( 2Lt0 *- | '° to ^^^. Thus the La- 
grangian case holds approximately when 

e 2L tat y > e (2Lt ot -U)y 1 



that is for 

2/max ~ In ( j- 1 ) (61) 

This estimate turns out to be quite good as we shall see. 

This result can be obtained also from the following argument. The largest 
y = |Im k\ that can be solution of Eq.(H) is approximately given by the equation 

1 - e^'-e^" '° i °A min = 

from where wc obtain the result of Eq. (pTj) . 

A similar argument was used by Kottos and Smilansky n& to obtain the 
gap empty of resonances < y < y m i n with 

1 1 ( l 

ymin ~ -, In 



'max \ | Amax | 



2.") 



In Subsection |6.2| , we presented a lower bound for this gap which is very 
accurate. Chaotic systems with a fractal set of trapped trajectories of partial 
Hausdorff dimension dn < 1/2, have a gap empty of resonances below the axis 
Re k given by 

y,„i„ = -P(l/2) 

This bound is based on the classical dynamics. In this case, the cumulative 
function h(y) vanishes for < y < y m i n = — P(l/2). 

The existence of the function h(y) in the limit x — Re fc — > oo for the case of 



graphs and the relations (|59|)-(|60|) are compatible with a conjecture by Sjostrand 
[£3J and Zworski |2J, g5| that the distribution of scattering resonances should 
obey a generalized Weyl law expressed in terms of the Minkowski dimension of 
the set of trapped trajectories, because this Minkowski dimension is equal to 
one for the quantum graphs. 

6.3.3 Width's distribution 

The density of resonances with a given imaginary part, P(y), is defined by 

{# of resonances k n = x n - iy n such that y < y n < y + dy} 

P(y)dy — - — r (62) 

{total # ol resonances} 

or in terms of the previously defined h(y): 

dh(y) 



P(y) 



dy 



The power law P(y) ~ y~ 3 ^ 2 has been conjectured in Ref. Q to be a 
generic feature of the density of resonances for systems with diffusive classical 
dynamics. The system studied in Ref. Q was the quantum kicked rotor whose 
classical limit is the standard map. For some values of the parameters that 
appear in this map, the phase space is filled with chaotic trajectories and not 
large quasi-periodic island are observed. Working with these particular values, 
the standard map produces a deterministic diffusion. In Ref. |lj, the kicked 
rotor has been turned into an open system by introducing absorbing boundary 
conditions at some fixed values, say +N/2 and —N/2. The distribution P(y) 
for the scattering resonances obtained in Ref. Jl2| is not identical to the one we 
observe for our graph. In particular, the density P(y) of Ref. 113] starts with 
P(y = 0) = and grow until a maximum value y ~ 7 c i, with j c \ the classical 
escape rate. It is the tail of P(y) which decreases from y w 7 c i as the power law 
P(y)~y~ 3/2 . 

The conjecture is that the power law holds for the tail of P(y), that is 
for large y. We present here the argument of Ref. (j^] which motivates this 
conjecture for systems with a diffusive classical limit. 

Consider a classical open diffusive system that extends from x = —N/2 to 
x = +N/2. Wc set absorbing boundary conditions at X = —N/2 and x = +N/2. 
Consider at t = a particle in the interval (—N/2,+N/2). Since the particle 
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escapes by a diffusion process the mean time that the particle takes to arrive at 
the border, starting at a distance X from it, is the diffusion time id ~ %■< We 
suppose that the resonant states are more or less uniformly distributed along the 
chain and that their quantum lifetime is proportional to the mean time taken by 
the particle to move from X to the border so that |Imfc| = Y ~ j- ~ -8%. The 
chain being symmetrical under x — > —x, we can assume a uniform distribution 
of X from the border x = —N/2 to the middle x = +N/2. The probability that 
the imaginary part is smaller that y is thus 



Prob {Y <y} = Prob 



aD 



< y \ = Prob 




<I = 1- 




where a is a dimensionless constant. Therefore, the probability density of the 
imaginary parts of the resonances k n — x n — iy n is given by 



P(y) = — Prob {Y<y} = 
ay 



N 



-y 



-3/2 



(63) 



This qualitative argument can be criticized on various points. In particular, 
the constant factor a is not determined and the assumption that the result holds 
for the quantum case is questionable because of the use of classical considera- 
tions. A complete theoretical validation of this law is thus lacking. However, 
the numerical results presented below give a support to this conjecture i n th e 
case of a multiconnected graph which is spatially extended (see Subsection 7.4 ) . 



7 Examples and discussion 

7.1 Simple graphs with two leads 

For the first two examples in this section, we consider Neumann boundary con- 
dition a 



ab 



In Figs. Ek and pic, we compare the decay of the quantum staying probability 
with the classical decay obtained from the Pollicott-Ruelle resonances and we 
also depict in Figs. Up and RH the spectrum of quantum scattering resonances 
for the corresponding graph. The initial wave packet also plotted in Figs, pp 
and [yd define a spectral window. 

The excellent agreement in Fig. |ja is due to the fact that the quantum and 
classical lifetimes of the resonances coincide. 

In fact, for the simple graph (inset Fig. [j]a), the scattering matrix is 



s = 


1 2 " 

2^ \ 
. 3 3 _ 


+ 9 + 3e 2ife s 


" 1 
1 


with lb = g. 


The resonances are given by the poles of S in Eq.( 


of 9 + Se 2lk 3 that is 




k n ± 


ln+1 In 3 

7T — I 

25 2 5 





(64) 



a), which are the zeros 



(65) 
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with n integer. These resonances have a lifetime r„ = p— with T n = — 2v„Im k n = 



,ln3 



where we used v = 2Re /c. This means that all the resonances have the 



lifetime given by 



win 3 



(66) 



This result could have been obtained using Eq.(p2[) and the classical Pollicott- 
Ruelle resonances from Eq. (f49|) . 
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Figure 1: (a) Decay of the quantum staying probability for the graph of the 
inset. The straight line is the exponential decay of the classical density for 
the same graph, (b) The dots represent the quantum scattering resonances of 
the graph in figure (a). The line superposed to the dots is the corresponding 
value associated with the classical decay rates of the classical Pollicott-Ruelle 
resonances. The curved line on the real axis is the initial wave packet, (c) The 
same as in (a) for the graph of the inset in (c). The steepest slope corresponds 
to the classical decay while the line with the smaller slope corresponds to the 
decay determined from the isolated resonance which is the closest to the real 
axis in figure (d). (d) The same as (b) for the graph of the inset in (c). 

In fact, for this graph, the only periodic orbit is the one that bounces on the 
bond b. The length of this periodic orbit is l p = 2g. The stability coefficient is 



2N 



a? a 1 - 

bb bb 



= 1 x | and we obtain X p l p — In 9. At 



given by \A p \ 2 = cxp(-X p l p ) 

the vertex 2 the particle is reflected with trivial back scattering and thus the 
analogue of the Maslov index fi p = 2 for the periodic orbit. Thus the quantum 
scattering resonances are the solutions of 

1 — exp I h 2igk + iix ] =0 



from which Eq. fl65|) follows, while the classical Pollicott-Ruelle resonances are 
given by the solutions of 

l-evp(bi9-2gs/v) = 

which follows from Eq.(f49|). These solutions are 

In 9 , 2inir .. 

•"-W^Wv (67) 

Therefore, all the Pollicott-Ruelle resonances have the lifetime t c \ = -p^ — 
—f—n- This lifetime coincides with the quantum lifetime obtained from the res- 
onances of the same graph in Eq. ( pq ) . 

The second graph in Figs. [Ip and |p has a chaotic classical dynamics. In 
this case, we also observe a very good agreement for short times between the 
quantum staying probability and the classical prediction and then a transition to 
a pure quantum regime. The decay in the quantum regime is again exponential 
because there is an isolated resonance that controls the long-time decay. In 
fact, as we see in Fig. |]d, there is an isolated resonance very near the real axis 
under the window given by the initial wave packet. The decay rate given by this 
resonance, T n — —4Rek n Imk n , gives the straight line with the smaller slope. 

7.2 Triangular graphs 

The following example is a nice illustration of the role played by the Lagrangian 



mean motion in the density of resonances described in Subsection 6.3. Consider 
a graph with the form of a triangle, that is we have three vertices and every 
vertex is connected to the two other vertices. The length of the bonds are a, b 
and c. See Fig. ||. 

Now, we add a semi-infinite scattering lead to each vertex and we use Neu- 
mann boundary condition at the vertices (see Fig. @a). In this case, the reso- 
nances are determined by the zeros of 

f(k) = e 2ifc ( a + b + c ) _ lQ e ik ( a + b + c ) _ g2ifc(a+6) _ e 2ik(b+c) 

_ e 2tk(a+c) _ 3 ( e ifca + jkb + gifcc) + 2 7 = 

replacing k = x — iy we have that the mean motion of f y {x) for y — > oo is 
M(y — » oo) = 2(a + b + c) and the density is given by h(oo) = a+b+c , which 
illustrates Eq. (159]). 
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(a) 



(b) 



Figure 2: The open triangle quantum graph: (a) One scattering lead is con- 
nected to each vertex, (b) Two scattering leads are connected to each vertex. 



Now, consider the same graph but with two semi-infinite leads attached to 
each vertex (see Fig. gb). In this case, the resonances are determined by the 
equation 



2ika 



I 2ikb I like , ik{a+b+c) 



replacing k = x — iy we have that the mean motion of f y (x) for y — > oo is given 
by max{a + b + c, 2a, 26, 2c} and therefore the density of resonances is given by 
h(oo) = maxj^r 2 , f , £, -f }, which illustrates Eq. @). 

An interesting observation is the following. The matrix T that contains the 
transmission and reflection amplitudes for the triangular graph can be written 
as 
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with u = — § — ' 3 and w = §■ + *■ ~ p; with /3 = 1 for the graph connected to 
two scattering leads at each vertex and (3 = when there is only one scattering 
lead per vertex. Note that det T = if /3 = 1 . We have computed the zeros 
of the function Zp(k) = det [I — T^D(fc)] for < (3 < 1. It is observed that 
some zeros in the lower half of the complex k plane decrease their imaginary 
part very fast as j3 increases. Therefore, we interpret the lowering of the density 
for the triangle connected to two scattering leads per vertex as the effect of 
resonances with |Imfc| — > oo. However, it should be noted that, even in the case 
that detT = 0, the analysis leading to Eq. ( pi] ) (i.e., the competition between 
the two leading terms of the function f v {x) for y — > oo) permits to obtain an 
approximate upper bound y max - 
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7.3 Spatially extended multiconnected graph 

Here, we consider the multiconnected graph of Fig. ||. The classical dynamics on 
this graph was studied in Ref. || where we showed that the escape is controlled 
by diffusion. 




Figure 3: Open graph built out of a periodic chain by attaching semi-infinite 
leads on a graph made of N unit cells. The figure shows a chain with N = 6 
unit cells. 



For this particular graph the T matrix is 



Tbb> = 



— I if the particle is reflected, i.e., b = b'; 

| for bonds b ^ b' which are connected; 

otherwise. 



The spectrum of quantum scattering resonances of this graph is depicted in Fig. 
for the chain with N = 5 unit cells. 
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Figure 4: Quantum scattering resonances for the graph in Fig. with N = 5 
unit cells. 



7.3.1 Width's distribution 

For this graph, we have computed the density P(y) of resonance widths defined 
by Eq. ( |62| ) . The histogram of resonance widths is plotted in Fig. |] for different 
sizes iV of the chain. Note that no resonance appears below y max = 2.276 which 
is the value computed from Eq.(|6l]). This value is independent of the system size 
N for N > 3. We have computed the eigenvalues of the matrix T for different 
values of N and we show the minimum eigenvalue in the following table: 
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Figure 5: P(y) in normal and log-log scales for the chain of Fig. || with 
(a),(t>) N = 6; (c),(d) N = 7. The continuous line shows the power law y~ 3 ^ 2 . 
The distribution stops at the value y max . The dashed line indicates the value 
y = |Imfc| = 7 c i/2 associated with the classical escape rate j c \. 
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I A m i n | 


yxnax 


1 


0.200 


4.553 


2 


0.392 


2.655 


3 


0.447 


2.276 



As we said for N > 3, this eigenvalue is independent of TV and is given by 
the case with N = 3. 

Figure [5] depicts the density of the imaginary parts in a log-log scale and the 
power law P{y) ~ y~ 3 ^ 2 conjectured in Ref. |12[ is observed. This is shown for 
chains of different lengths in these figures. The power law holds for not to small 
sizes N. This observation is therefore a support for the conjecture of Ref. ]12j 
in the case of multiconnected graphs as the one of Fig. H. 
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7.3.2 Detailed structures in the resonance spectrum 

Let us now discuss the details of the distribution of resonances. 

First of all, we notice that the statistical description we have considered, that 
is the distribution of y, is based on the homogeneous character of the resonance 
distribution along the x — Re k axis. This homogeneity indeed holds on large 
scales as we see in Fig. || Nevertheless, at small scales, we can see in Fig. g the 
formation of bands of quantum resonances characteristic of a periodic system 
for which Bloch theorem applies pa. The reason is that, for long times, the 
system is able to resolve finer scales in wavenumber (or energy) and thus after a 
given time the system will "feel" the periodic structure and evolves ballistically 
as follows from Bloch theorem. 



S 0.7 



Re it 

Figure 6: Quantum scattering resonances (dots) and transmission probability 
amplitude T sgsi (solid line) for the chain of Fig. || with N — 7 unit cells. 

The distribution P(y) shows a peak at the origin (i.e. y ~ 0) which grows 
with the number of unit cells, N. This behavior tells us that "bands" of reso- 
nances with small y are created when we increase the size of the chain. In fact 
we have that, as for the resonances of one-dimensional open periodic potentials 
[Em, we have TV — 1 resonances per band and that the bands converge as ~ 1/N 
to the real axis when N increases. This is consistent with the ballistic behavior 
that should be observed in the long-time limit. 

On the other hand, the resonances with larger values of y = |Imfc| are 
not arranged in a band structure and their number does not increase when we 
change the size N of the system. In fact, these resonances are located at the 
same position in the complex plane for every value of N. Therefore we interpret 
such short-lived resonances as metastable states that decay without exploring 
the whole system. 

In Fig. raa, we superpose the resonance spectrum for chains with two different 
sizes. The figure shows that, for resonances with large values of y — |Imfc|, 
neither their position, nor their number change with N, but resonances with 
small values of |Im/c| converges to the real axis as we increase N. Therefore, 
the relative number of resonances in the tail of the distribution decreases as 
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compared to the number of resonances with small values of y which increases 
with TV. We can conclude that the resonance spectrum converges to the real 
axis in probability when we increase N. This behavior is in contrast with the 
simple systems analyzed in Ref. p™ , for which we observed that each resonance 
converges to the real axis as N increases. 

Eq.(|63|) predicts a relative decrease of P(y) as 1/N. This law is verified as 
shown in Fig. [fl where we plot for clarity the function N x P(Im k) for only a 
few cases. According to the theoretical distribution of Eq.(|63]) N x P(Imfc) oc 
V^DlImfcl -3 / 2 , where the right-hand side is independent of N and moreover, 
the proportionality factor is determined by the diffusion coefficient of the chain. 
We have computed the diffusion coefficient for this graph in Ref. ||. The 
continuous straight line in Fig. corresponds to the diffusion coefficient of the 
infinite chain. The good agreement shows that the proportionality constant a 
in Eq. (63) is of order one, as expected. 




Figure 7: Superposition of the log-log plots of N x P(|Imfc|) as a function of 
y = |Im fc| for N = 5, 6, 8 and the same graph as in Fig. 0. 

In Fig. pi we observe deviations from the power law at small and large values 
of |Imfc|. The dashed lines in Fig. a indicate the value corresponding to the 
classical escape rate. We see from the figure that the distribution of imaginary 
parts is well described by the power law for y = |Imfc| > 7ci/2 and the classical 
escape rate 7 c i is near the transition to this power law. On the other hand, 
the tail of the distribution also deviates from the power law. The distribution 
follows the power law until a value which decreases when N increases. Beyond 
this value the distribution seems to fluctuate around an almost constant value, 
and then drops rapidly to zero. 

The region where the distribution of imaginary parts fluctuates around some 
value corresponds to the region where the resonances are independent (in num- 
ber and in position) of the value of N. This can be seen in Fig. || In Fig. ^a, 
we depict the resonances for N — 7 and N = 8 and we draw a line which sepa- 
rate resonances that belongs to bands and change with N from those that are 
independent of N. This line is given by y = |ImA;| = 0.5. In Fig. pp, we depict 
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Figure 8: (a) Resonances for the chain of Fig. |^ with N = 7 (circles) and N = 8 
(crosses), (b) The distribution of resonance widths for the chain with N = 7. 
The straight line in figure (a) and (b) separates the two families of resonances: 
those which depend on N from those which do not. This figure shows that the 
change in the distribution of the imaginary parts of the resonances is associated 
with the different families of resonances. 

the density of resonances for N = 7, the value y = |Imfe| = 0.5 is indicated by 
the vertical straight line. This line marks the separation between both families 
of resonances and thus the limit of the power law y~ 3 ' 2 . 

7.4 Linear graph with emerging diffusion 

For extended periodic open graphs, the classical decay given by the leading 
Pollicott-Ruelle resonance, corresponds to the decay of a diffusion process M. 
Here, we analyze the decay of the quantum staying probability for such a graph. 
Here, we consider a linear periodic graph with a unit cell composed by two 
bonds a and b of incommensurate lengths l a and If,, respectively. At the vertex 
that join these two bonds we have a scattering matrix 17(772) and the vertices that 
join two unit cells have a scattering matrix <j(?7i). These scattering matrices are 
of the form 



<j{rj) 



1 sm 77 cos rj 
cos rj i sin rj 



(68) 



Therefore, the transmission and reflection probabilities for the classical dynam- 
ics are Ti = cos 2 (r/i) and Ri — sin 2 (77,). 

Fig. depicts an open graph by considering only N units cells connected to 
semi-infinite leads at the left-hand and right-hand side of the finite graph. 

7.4.1 Classical diffusive behavior 

In the limit N — > 00 of an infinitely extended chain, the graph becomes periodic 
and the motion is diffusive. The diffusion coefficient can be calculated by con- 
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Figure 9: Open graph built out of a periodic linear chain with a unit cell made 
of two bonds with scattering matrices (j^) at each vertex. The figure shows a 
chain with N = 2 unit cells. 

sidering the Frobenius-Perron-type matrix Q(s, q) defined in Ref. M in terms 
of the classical wavenumber q and the rate s. For the present graph, we have 
that 



Q(*,g) 






e +«2 e -s^Ti 




e" s ^i?i 




e- s ^T 2 










e 


e- s ^R 2 










e 





e « Hi 


e 


-ii e -^Ti 







« il 2 

S ^T 2 




(69) 



where the columns and rows are arranged in the following order (a, b, a, 6). The 
diffusion coefficient is obtained from the second derivative of the first branch at 
q = 0. Developing det[l — Q(s, q)] for small values of q and s we get 



det[l - Q(s, q)] = 2/(i?!T 2 + Tii? 2 )- + TiT 2(7 2 + 0(s 2 ) + 0{sq 2 ) 

v 



0(q 4 



The diffusion coefficient D is thus given by 



D _ i gVg) 

2 9 9 2 



uTiT 2 



q=0 



2/(i?iT 2 +Tii? 2 ) 



(70) 



where / = / a + If, is the total length of the unit cell. 

In Ref. H , we have considered other examples and we have shown that the 
classical escape rate 7 c i(AT) for a finite open chain of size N is well approximated 
by Eq. (51) in the limit N — > oo. For the present linear graph, Fig. [l(] illustrates 



that the classical lifetimes of the chain are indeed determined by the diffusion 
coefficient according to t c \ = l/7ci ~ N 2 in the limit N — > oo. 



7.4.2 Spectrum of scattering resonances and its gap 

The resonance spectrum is depicted in Fig. [Ty for different chain sizes N = 1, 
2, 3, and 4. 

The structures of the resonance spectrum and the presence of a gap for the 
sizes N = 1 and N — 2 can be understood thanks to the topological pressure 
plotted in Fig. [12 for these graphs. We see that the chains with N = 1 and N — 
2 unit cells has a value P(l/2) < and, therefore, a gap empty of resonances. 
This gap appears in Figs. Oa and Ob. For the chains with N = 3 and N = 4 
unit cells, P(l/2) > and we do not have this upper bound for Figs, flp and 
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Figure 10: Plot of the lifetime r c i(iV) = l/~/ c i(N) for the chain of Fig. | with 
AT = 1, 2, 3 and 4 unit cells as a function of TV 2 . The dashed line corresponding 
to the diffusion coefficient D — 0.506i> given by Eq.(70) is approached by the 
classical lifetimes for increasing size N of the open graphs. We have used the 
parameters 771 = 0.1 and r\i = (\/5 — l)/2 and l a = 0.5 and lb = \[2. 



7.4.3 Decay of the quantum staying probability 

For these open graphs, we have calculated the quantum staying probability in 
order to illustrate the emergence of a diffusive behavior. 

The quantum time evolution has been computed by using the method de- 
scribed in Subsection 5.4. The propagator was calculated by fast Fourier trans- 
form from the Green function. The staying probability was calculated by con- 
sidering an initial Gaussian wave packet located on a bond 1 inside the graph: 



My,o) = 



1 



7T cr(l - e 



-2E<7 2 



1/2 



sm(v£,y)e i"' 2 



This wave packet is built by the modulation of a Gaussian with a sine function 
of wavenumber v E. The wave packet is thus centered at y = yo (and we take 
2/0 = ii/2). The wave packet is effectively on the bond 1 if Ay = 2rr « l\. 

In Fig. [13L we compare the decay of the quantum staying probability with 
the decay obtained from the leading Pollicott-Ruelle resonance for different sizes 
N (see Fig. [Tc]). This leading Pollicott-Ruelle resonance is the classical escape 
rate sq = — 7cl(^) at the velocity corresponding to the mean energy E of the 
Gaussian wave packet. The agreement observed in Fig. O for short times 
between the classical decay and the quantum decay shows that for short times, 
the quantum evolution follows a diffusion process. 
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The deviation that appears in Fig. O for longer times corresponds again 
to the decay determined from an isolated resonance and is therefore a pure 
quantum effect. 

8 Conclusion 

We have studied dynamical and spectral properties of open quantum graphs 
with emphasis in their relation to the transport properties of the corresponding 
classical dynamics. The classical dynamics has been obtained as the classical 
limit of the quantum dynamics. 

The time evolution of the quantum system is obtained from the propagator 
which is computed as the Fourier transform of the Green function. A closed 
expression and a multi-scattering representation for the Green function has been 
presented. We want to emphasize that both, classical and quantum evolution are 
considered with the continuous time in opposition to the discrete time evolution 
often considered in the literature of lattice networks. 

In particular, we have computed the quantum staying probability. For short 
times, this quantity decay exponentially in time with the classical escape rate. 
This classical escape rate is obtained from the classical zeta function and corre- 
spond to the leading Pollicott-Ruelle resonance. For large open periodic graphs 
the leading Pollicott-Ruelle resonance determines a decay dominated by diffu- 
sion and the decay of the quantum staying probability reveals at the quantum 
level the emerging diffusion process. 

On the other hand, the quantum spectral properties are also related to trans- 
port. The resonance spectrum reveals some features related to both, the ballistic 
propagation in the periodic system for the long time evolution and the diffusive 
classical dynamics that is an approximation for the quantum dynamics for times 
shorter than the Heisenberg times. 

Indeed, the Fourier transform relates the long-time behavior to the variations 
at small energy scales. At small scales, the resonance spectrum is arranged in 
bands of N — 1 resonances that converges as ~ 1/N towards the real axis. 
Therefore the lifetime of each resonance in a band is proportional to the size of 
the system reflecting the ballistic transport that characterize periodic system as 
state the Bloch theorem [g6| . 

Nevertheless, this ballistic propagation affects the long-time dynamics. In- 
deed, the Bloch theorem gives information about stationary states. The short 
time dynamics is determined by the large fluctuations of the resonance spec- 
trum. At a large scale, the resonance spectrum is uniform over the Re k axis 
(see Fig. I|). This distribution display a power law P(y) ~ y~ 3 / 2 early obtained 
in Ref. | 12[ and conjectured to be a general law for open quantum systems 
with diffusive classical limit. Our numerical results support this conjecture for 
multiconnected spatially extended graphs. 

Moreover, we have presented an alternative derivation of the density of reso- 
nances obtained in Ref. U& , which is based on general properties of the secular 
equations, namely, that is an almost-periodic function which in the appropriate 
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limit its mean motion is in the Lagrangian case. This method allows us to ob- 
tain an approximate lower bound for the distribution of y = |Imfc|. Moroever, 
an upper bound for this distribution has been obtained from the topological 
pressure P(/3). This bound creates a gap in the resonance spectrum under the 
condition P(l/2) < and is absent otherwise. 

In summary, we have studied quantum properties of extended open periodic 
graphs. We have shown that the Pollicott-Ruelle resonances determine the 
decay of the quantum staying probability, which in turn shows the appearance 
of diffusion at the quantum level. The diffusion process is thus encoded in the 
distribution of scattering resonances. 
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Figure 11: Scattering resonances of the chain of Fig. O with N = 1, 2, 3 and 
4 unit cells. We have used the parameters i]i = 0.1 and r/ 2 = (\/5 — l)/2 and 
l a = 0.5 and lb = V%- 
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Figure 12: From bottom to top, we plotted the topological pressure for the chain 
of Fig. g with the same parameter values as in Fig. [ll] and the sizes N = 1, 2, 
3 and 4 unit cells, respectively. 
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Figure 13: From bottom to top: Decay of the quantum staying probability for 
the chain of Fig. [| with 1,2,3 and 4 unit cells and the same parameter values as 
in Fig. [jl]. The straight lines gives the corresponding classical decay as obtained 
from the Pollicott-Ruclle resonances. 
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